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Abstract. A new class of solutions of the Einstein field equations in spherical 
symmetry is found. The new solutions are mathematically described as the met- 
' rics admitting separation of variables in area-radius coordinates. Physically, they 

, describe the gravitational collapse of a class of anisotropic elastic materials. Stan- 

dard requirements of physical acceptability are satisfied, in particular, existence 
of an equation of state in closed form, weak energy condition, and existence of 
a regular Cauchy surface at which the collapse begins. The matter properties are 
generic in the sense that both the radial and the tangential stresses are non van- 
ishing, and the kinematical properties are generic as well, since shear, expansion, 
I and acceleration are also non-vanishing. As a test-bed for cosmic censorship, 

the nature of the future singularity forming at the center is analyzed as an exis- 
tence problem for o.d.e. at a singular point using techniques based on comparison 
(•~^ ' theorems, and the spectrum of endstates - blackholes or naked singularities - is 

CO ! found in full generality. Consequences of these results on the Cosmic Censorship 

' conjecture are discussed. 

o 

O 

1. Introduction 

! The final state of gravitational collapse is an important open issue of classical 

gravity. It is, in fact, commonly believed that a collapsing star that it is unable to 
radiate away - via e.g. supernova explosion - a sufficient amount of mass to fall 
^ . below the neutron star limit, will certainly and inevitably form a black hole, so 

that the singularity corresponding to diverging values of energy and stresses will 
be safely hidden - at least to faraway observers - by an event horizon. However, 
this is nothing more than a conjecture - what Roger Penrose first called a "Cosmic 



Censorship" conjecture - and has never been proved. On the contrary, in the 
last twenty years of research, many analytic examples of (spherically symmetric) 
naked singularities satisfying the principles of physical reasonableness have been 
discovered. 

Spherically symmetric naked singularities can be divided into two groups: those 
occurring in scalar fields models ^ and those occurring in astrophysical sources 
modeled with continuous media, which are of exclusive interest here (see [ |TT| ] for a 
recent review). The first (shell focusing) examples of naked singularities where dis- 
covered numerically by Eardley and Smarr [Q] and, analytically, by Christodoulou 
[01, in the gravitational collapse of dust clouds. Since then, the dust models have 
been developed in full details: today we know the complete spectrum of endstates 
of the gravitational collapse of spherical dust with arbitrary initial data [|T^. This 
spectrum can be described as follows: given the initial density and velocity of the 
dust, a integer n can be introduced, in such a way that if n equals one or two, the 
singularity is naked, if n > 4 the singularity is covered, while if n = 3 the system 
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undergoes a transition form naked singularities to blackholes in dependence of the 
value of a certain parameter. If the collapse is marginally bound, the integer n is 
the order of the first non vanishing derivative of the energy density at the center. 

The dust models can, of course, be strongly criticized from the physical point of 
view. Since from one side very few results are known on perfect (i.e. isotropic) flu- 
ids and, on the other side, stresses have to be expected to be anisotropic in strongly 
collapsed objects, one of us initiated some years ago a program whose objective 
was to understand what happens if the gravitational collapse occurs in presence of 
only one non- vanishing stress, the tangential one [|2^, |21|]. The program can be 
said to have been completed: we know the complete spectrum of the gravitational 
collapse with tangential stresses in dependence of the data both for the clusters 
[ |T0| , and for the continuous media models 

The results of this program are somewhat puzzling. In fact, what one should 
expect on physical grounds is that the equation of state plays a relevant role in 
deciding the final state of the collapse. Well, it is only apparently so. In fact, also 
in the tangential stress case a parameter n can be defined, in such a way that the 
endstates of collapse depend on n exactly in the same way as in dust. 

In the present paper we present a complete, new model of gravitational collapse 
which includes both radial and tangential stresses. This is done deriving a class of 
anisotropic solutions which is in itself new, and contains the dust and the tangential 
stress metrics as special cases. As far as we are aware, this is the first time that the 
spectrum of the endstates of a solution satisfying all the requirements of physical 
reasonability and exhibiting both radial and tangential stresses is found. 

As will be discussed in the concluding section, in view of the results of the 
present paper, the case for a cosmic censor - at least in spherical symmetry - be- 
comes very weak. 



2. Einstein Field Equations in area-radius coordinates 

Consider a spherically symmetric collapsing object. The general line element in 
comoving coordinates can be written as 

(2.1) ds^ = -e^'^dt^ + (l/r^) dr^ + R^^dO'^ + sin^ ^d^^) 

(where z/, rj and R are functions of r and t). We shall use a dot and a prime to denote 
derivatives with respect to t and r respectively. 

In the present paper we consider as admissible sources of the gravitational field 
only those matter models which admit a well-defined thermodynamical description 
in terms of the standard relativistic mechanics of continua. We recall that the phys- 
ical properties of an elastic material in isothermal conditions can be described in 
terms of one state function w, which depends on three parameters (the so called 
strain parameters) and on the coordinates. In the comoving description, the strain 
tensor of the material can be expressed in a "purely gravitational" way and, as a 
consequence, the state function depends only on the space-space part of the met- 
ric(see e.g. [JIP]] and references therein). In addition, in spherical symmetry, the 



state function cannot depend on angles, and therefore the equation of state of a 
general spherically symmetric material can be given as a function of r and of two 
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Strain parameters : 
(2.2) 



w = w{r, R, T]) 



It is useful also to introduce the matter density 
(2.3) p- ^ 



(where E = E{r) is an arbitrary positive function) so that the internal energy 
density is given hy e = pw. If w depends only on p, one recovers the case of the 
barotropic perfect fluid; in general, however, the stresses are anisotropic and are 
given by the following stress- strain relations: 



(2.4) 



dw 
dr] ' 



Pt 



1 ^dw 



Comoving coordinates r, t are extremely useful in dealing with gravitational col- 
lapse because of the transparent physical meaning of the comoving time. We shall, 
however, make systematic use here also of another system of coordinates, the so- 
called area-radius coordinates, which were first introduced by Ori p3| ] to study 
charged dust, and then successfully applied to other models of gravitational col- 
lapse (see e.g. [jT^, ^). These coordinates prove extremely useful for technical 
purposes, as will be clear below. In these coordinates the comoving time is replaced 
by R. The velocity field of the material = e~^5f transforms to v'^ = e~'^RS^ 
and therefore the transformed metric, although non diagonal, is still comoving. It 
can be written as 

(2.5) ds^ = -Adr^ - 2BdRdr - CdR^ + R^{de^ + sin^Odip'^) , 
where A, B and C are functions of r and R. Obviously 



(2.6) 



where we have denoted u = \Re~'^\. Formula ( |2.2D implies that the internal energy 
w now depends on two coordinates (r and R) and on only one field variable r]. It is 
convenient to introduce the quantity 

1 



(2.7) 



A :=5^ 



AC 



> 



(so that y/r] = 1/{u^/A) and B = — a/A + A/u'^), and to use A,u and as 
the fundamental field variables. A convenient set of Einstein equations for these is 
G; = SttT;, = SttT/^ and = 8nT^ 0. Denoting partial derivatives with 
a comma, the first of these equations reads 



1 

R^ 



(2.8) 

so that A decouples from u and A: 
(2.9) A = A 



A [A 
A V A 



2E(r) 2A 



R 



R 
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where 



(2.10) 



A : = 



R 



A'Ka'^Pr{r, 0", //(r, a)) da. 



Ro{r) 



In the above formulae, the arbitrary positive function RQ{r) describes the values of 
R on the initial data and will be conveniently taken to be 

(2.11) i?o(r) = r. 

To find the physical meaning of the other arbitrary function F, we introduce the 
so called Misner-Sharp mass \l'(r, R). It is defined in such a way that the equation 
R = 2^ spans the boundary of the trapped region, i.e. the region in which outgoing 
null rays reconverge. This boundary is the apparent horizon of the spacetime, and 
it can be shown that (see e.g. 0]): 

(2.12) ^(r,i?) = {R/2) (1 - gf^'d^Rd^R) . 
Transforming to radius-area coordinates and using ([2.9[), one easily gets 

(2.13) ^ = (i?/2) (1 - A/A) = F(r) + A 

and therefore the function F is the value of the Misner-Sharp mass on the data. 
Since A vanishes if the radial stress is zero, when = the Misner-Sharp mass 
is conserved during the evolution (i.e. independent on R). This is the main techni- 
cal reason which makes the spacetimes with vanishing radial stresses simpler with 
respect to general collapse models. 

Using (Pj^) it can be shown that the two remaining field equations can be written 
as follows: 



(2.14) 

and 

(2.15) 



orj J 



2m 



2* 
R 



U 



3. The general solution admitting separation of variables 
3.1. The solution. Let us consider the system of two coupled PDE's for A and u 



( |2.14[ )-( p.l5D . If A does not depend on the field variables, equation ( |2.14{ ) becomes 
algebraic and the system decouples. Functional properties of A are related to the 
choice of w and in fact, A will be independent of rj if pr does. This obviously hap- 
pens if Pr is zero (leading to the already well known cases of dust and of vanishing 
radial stresses, for which w does not depend on rf) but also if p^ is some function of 
r and R only. This in turn occurs if w is of the form 



(3.1) 



w{r,R,ri) = h{r,R) + —l{r,R), 



the case Pr = corresponding to / 
non-vanishing 



0. If / is non zero the radial pressure is 
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SO that 
(3.3) 

Using (jXTj) and ( Pl] ), equation ( P^ gives 



A = / /(r, a) da. 



2^ 

E^i/,r = h\lu^ + l-—, 



which allows us to compute the quantity u^: 

(3.4) u^{r,R) = 



2^ fE'^r 



As a consequence, ( |2.15p becomes a quadrature which allows the calculation of 



(3.5) VA 
where we have defined 

(3.6) n{r, R) = 



1 dn ^ 

Y(r, a) dr 



r, a) da, 



R 



2^(r, R) + R{Y^{r, R) - I)' 



(3.7) 



r(r, R) 



E(r)^,^(r, R) 
h{r, R) 



It is convenient to eliminate the indefinite integral in (p.Sp. For this aim, it is useful 
the following equation 



(3.8) 



Y{r,R) = R'ir,tir,R))^{r,tir,R))., 



that can be easily found from the relations between the metrics coefficients in the 
comoving and the area-radius coordinates. Using (p77|), ( p.8D and the fact that 



R' = Ifor R = r, we get \/A(r, r) = ^7^, which fixes the integration constant, 
to obtain 

0.9) ^{r, R) = /" y^.^ir, -)^- + ^- 

J n y (r, cr) or r (r, r) 

3.2. Physical requirements. In this section we will discuss the conditions that 
have to be imposed on the equation of state and on the data, in order for the solu- 
tions to be physically meaningful. 

From now on, we consider only spacetimes whose matter source satisfies the 
equation of state (|3.1D. For such spacetimes, using ( |2.13| ) and (|3.3D, the constitutive 
function w may be written as 



(3.10) 



w{r, R, 7]) = E[r) + — 

V Y[r, R) ^ 
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therefore the energy density e = pw and the stresses (E^ are given by: 



(3-11) -4^(#^^^ + ^^- 



(3.12) Pr 



(3-13) = -77^ -TT T7VtT7T + 



87ri? \ Y OR ^ 

From now on we choose as independent arbitrary functions \& and Y. In the 
present paper, we assume C'^ regularity of the data: the equation of state and the 
arbitrary functions are assumed to be Taylor-expandable up to the required order. 
Due to eq. ( |3.1 ID , ( |3.12D and ( [3.1 3D , it will then be possible to translate all the 



conditions to which the metrics have to satisfy in order to be physically meaningful 
(like e.g. weak energy condition or regularity of the center up to singularity for- 
mation, see below), in terms of conditions on the functions ^ and Y or on their 
derivatives. 

We first impose the weak energy condition (w.e.c). In spherical symmetry, this 
condition is equivalent to e > 0, e + pr > 0, e + pt > 0. These inequalities lead to 

(3.14) v]/, >o, VI/ ^>o, VI/, >|r(^^^ , ^^n>^^f^RR. 

For an explicit solution satisfying weak energy condition we refer to the example 
in ^ below. 

We also impose regularity of the metric at the center ('local flatness'). In comov- 
ing coordinates r, t this amounts to require 

(3.15) R{0, t) = 0, e^(°'*) = i?'(0, t). 

In addition, the stress tensor must be isotropic at r = 0, that is 

(3.16) pr{0,t)=ptiO,t), 

for any regular t = const hypersurface. Finally, we require the existence of a 
regular Cauchy surface (t = 0, say) carrying the initial data for the fields. These 
requirements are fundamental, since they assure that the singularities eventually 
forming will be a genuine outcome of the dynamics. 

We have already chosen R = r at t = 0; using ( P.llD , ( p^ and the fact that 
R'{r, 0) = 1 we find the expression for the initial energy density eo(r): 

(3.17) ^__|^|^4V(r,r) + *^«(.r,r) 



47rr^ 

Once the regularity conditions are satisfied, the data will be regular if this func- 
tion is regular. For physical reasonability we also require the initial density to be 
decreasing outwards, that is eQ(r) < for r > 0. It is easy to check, that the 
above stated conditions are equivalent, in terms of the functions ^ and Y, to the 
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following: 
(3.18) 

F(0,0) = 1, 
(3.19) 

^(0,0) = *,,(0,0) = ^,fi(0,0) 
(3.20) 



^,rr(0,0) 



^,,,j(0,0)=^,fi«(0,0) = 0, 



^,rr{r, r) + 2^,rRir, r) + ^Rnir, r) - - (^,,(r, r) + ^,^(r, r)) < 0. 

r 

Definition 3.1. We say that a spacetime is a physically valid area-radius separable 
spacetime (ARS) if 

(1) the equation of state of the matter is ( |3.1(J| ), where \l/(r, R) and Y{r, R) are 
arbitrary positive functions; 

(2) the weak energy condition, the regularity condition, and the condition of 
decreasing initial density hold, i.e. ^{r, R) and Y{r, R) satisfy ( p.l4t ) and 

3.3. Kinematics. Spherically symmetric non-static solutions can be invariantly 
classified in terms of their kinematical properties (see e.g. [[IT], [TEp). Most of 
the known solutions have vanishing shear, or expansion, or acceleration. For the 
solutions studied here all such parameters are generally non-vanishing. 

The acceleration in comoving coordinates r, t is given by = z/'(5^, therefore, 
it can be uniquely characterized by the scalar A := ^a^a^. Using ( p?8l ) together 
with the relation 

(3.21) Yv' = R'Yr 

(which easily follows from the the field equation in comoving coordinates R' = 
Rv' + R'X) we get 

(3.22) A = Y^R. 
The expansion 6 is given by 



(3.23) 



Q = ±u 



u^fK R 



Finally, the shear tensor cr^ can be uniquely characterized by the scalar a 



(2/3) ^a^^a- given by 
(3.24) 



3 



1 

R 



3.4. Special classes. Special cases of the solutions discussed above are: 

(1) The dust (Tolman-Bondi) spacetimes. The energy density equals the matter 
density, and this implies = and Y = -E\l/,r- 

(2) The general solution with vanishing radial stresses. Vanishing of pr implies 
^ .R = (see ( |3.12| )), while Y depends also on R (if Yr = one recovers 
dust). The properties of these solutions have been widely discussed in [^0]]. 
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(3) An interesting new subclass is obtained by imposing the vanishing of the 
acceleration. This subclass corresponds to a choice of the function Y de- 
pending on r only. It is well known that acceleration-free perfect fluid mod- 
els can describe only very special collapsing objects (the pressure must be 
a function of the comoving time only) while these anisotropic solutions ex- 
hibit several interesting features, like e.g. a complete spectrum of endstates 

3.5. An explicit example: "anisotropysations" of Tolman-Bondi-de Sitter spa- 
cetime. The investigation of various explicit examples of new solutions within the 
class presented here, as well as their physical applications, will be presented in a 
forthcoming paper [0]. We restrict ourself here only to sketch a simple acceleration 
free (i.e. Y = Y(r)) case, essentially with the aim of showing explicitly that the 
new sector of the solutions (w.r. to dust and tangential stress case) is non empty. 

A simple way to produce explicit examples of new collapsing solutions within 
the equation of state considered here is to choose the mass function as the sum of 
a function of r only (governing the dust limit) and a function of R only (governing 
the anisotropic stresses). In order to satisfy the requirements of Definition pTT] , easy 
calculations show that the function ^ must have the form 

(3.25) ^(r,i?)= / 0(s)s^ds+ / xM^^da, 

Jo Jo 

where x ^^id are positive and non increasing functions. In the particular case 
in which x is a constant, the contribution of x to the energy-momentum tensor is 
formally the same as that of a cosmological constant, and therefore the spacetime 
coincides with the so called Tolman-Bondi-de Sitter (TBdS), describing the col- 
lapse of spherical dust. This is actually the unique model of gravitational collapse 
in presence of a lambda term which is known so far. In recent years, the descrip- 
tion of gravitational collapse in such spacetimes attracted a renewed interest both 
from the astrophysical point of view, since recent observations of high-redshift type 
la supemovae suggest a non-vanishing value of lambda, and from the theoretical 
point of view, after the proposal of the so-called Ads-Cft correspondence in string 



theory. In such a context, the solutions ( |3.25D can be used to investigate the effects 



of stresses by adding higher order terms to the function x- For instance, choosing 

v[/(r, i?) = ar^ + 



1 + R/R, 

where R^ is a positive constant, one obtains a model which is TBdS homogeneous 
and isotropic near the center but becomes anisotropic and inhomogeneous when R 
is of the order of R^ (we stress, however, that this is only one possible application 
of this special sub-class of solutions). 

4. Conditions for singularity formation 



4.1. Shell crossing and shell-focusing singularities. Due to eqs.(p7l) and ( [3.1 ID , 



the energy density becomes singular if, during the evolution, R or uV A vanish. 
The latter case corresponds to R'{r, t) = in comoving coordinates and is called a 
shell-crossing singularity, since it is generated by shells of matter intersecting each 
other. Shell-crossing singularities correspond to weak (in Tipler sense) divergences 
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of the invariants of the Riemann tensor. As a consequence, these singularities are 
usually considered as "not interesting", although no proof of extensibility is as yet 
present in the literature. In any case, we shall concentrate here on the the shell- 
focusing R = singularity, for which no kind of extension is possible. Therefore, 
we work out the conditions ensuring that crossing singularities do not occur. 

Since we are considering a collapsing scenario, R < r during evolution. There- 
fore, due to eq. (|J!9|), it is sufficient to require strict positivity of the function 
-\/A(r, 0) together with non-increasing behavior of \/A{r, R) w.r. to R. We thus 
have the following: 

Proposition 4.1. In a ARS spacetime the formation of shell-crossing singularities 
does not occur if 

(4.1) m^r{r,R) + Y{r,R)Y^r{r,R)R>{i /or r > 0, i? G [0, r], 



(4.2) rw^?^(^'^)^^+W^>0 /orr>0. 
Jo y[r,(j) or Y{r,r) 

The above condition, although being only sufficient, allows for a wide class of 
new metrics. It becomes also necessary in the particular case of dust [[131]. 

Spherically symmetric matter models do not, of course, unavoidably form sin- 
gularities if stresses are present. One can, in fact, construct models of oscillating or 
bouncing spheres. In comoving coordinates (r, t), the locus of the zeroes of R{r, t) 
- if any - defines implicitly a singularity curve ts{r) via R{r,ts{r)) = 0. The 
quantity ts(r) represents the comoving time at which the shell labelled r becomes 
singular, and therefore the central singularity forms if 

(4.3) lim tsir) = t,, 

is finite (t* is positive due to the regularity of the data at t = ). To study the 
behavior of this limit in dependence of the choice of the data, we use R = —e^u. 
Integrating along a flow line we get 

(4.4) t= / e-^("'")H(r, a)da, 

where the initial condition ( [2.1 1| ) has been used. Up to time reparameterizations we 
can assume z/(0, t) = 0. Using this fact, together with ( [3.21[ ), it can be seen that u 



is uniformly bounded and limr^o+ '^(^i cr) = uniformly for a G [0, r]. Therefore 
it is necessary and sufficient to require that, for r near 0, 7Y(r, a) is integrable with 
respect to a in [0, r] and 



(4.5) < lim / n{r, a)da < oo. 

'"^0+ Jo 

For this aim, let us consider the following Taylor expansion centered at the point 



(0, 0), where we take Definition |3TT[ into account: 
(4.6) 

2'i'{r,R) + R{Y\r,R) -1) = 2Yr{0,0)rR + 2YR{0,0)R^ + J2 Yl K^'^^ 

k>3 i+j=k 
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Changing variable a = rr and recalling the definition of H (eq. (|3.6D) the integral 
in ( |4.5p becomes 

= dr, 

2r,(0, 0)r + 2Yr{0, 0)r2 + r ^.+,=3 hijT^ + ^ v{r) 

with Lp(0) = 0. Convergence of this integral to a finite non-zero value requires 
vanishing of the zero order terms and at least one non vanishing first order term in 
the denominator (i.e. {h^o, /i2i, hi2) 7^ (0, 0, 0)). We shall, however, consider only 
the sub-case in which h-^o (which we will call a hereafter) is non-vanishing, since 
a vanishing a would correspond to a bad-behaved "dust limit" of the solutions, i.e. 
when the source of the anisotropic stress becomes very weak. For the same reason, 
a cannot be negative otherwise the weak energy condition would be violated in the 
same limit. Thus: 

Proposition 4.2. In an ARS spacetime shell focusing singularities form if 
(4.7) F,(0,0) = y^(0,0) = 0, a>0. 

From the above discussion, we now define the subclass of ARS spacetimes that 
we are dealing with hereafter. 

Definition 4.3. An ARS spacetime is called collapsing if: 

(1) Shell-crossing singularities do not form (i.e. ^ and Y satisfy to (|4.1D, (|4.2D); 



(2) Shell-focusing singularities form in a finite amount of comoving time (i.e. 



and Y satisfy to (p^D). 



4.2. The apparent horizon and the nature of the possible singularities. A key 

role in the study of the nature of a singularity is played by the apparent horizon (see 
for instance [fT^]). The apparent horizon in comoving coordinates is the curve th{r) 
defined by R{r,th{r)) = 2\E'(r, i?(r, t/i(r))). In area-radius coordinates, one has 
correspondingly a curve Rh{r) and, since \E',i?(0, 0) = 0, implicit function theorem 
ensures that Rh{r) is well defined in a right neighborhood of r = and such that 
Rhir) < r. This fact, along with the requirement ^ « > coming from the w.e.c. 



( |3.14| ) ensures \E'(r, Rhir)) < ^(r, r). Now using ( |4.4| ) it is 



(4.8) < ts{r) - th{r) = e-^Hir, a)da < e-'"Hir,a)da 

Jo Jo 

Changing variable a = rr as before and using the fact that u is bounded and 
\E'(0, 0) = 0, we find that this integral converges to as r tends to 0, ensuring the 
following 

Proposition 4.4. In a collapsing ARS spacetime the center becomes trapped at the 
same comoving time at which it becomes singular, that is 

(4.9) lim th{r) = U= lim ts{r). 

r^0+ r-»0+ 

In the models studied in this paper, the only singularity that can be naked is the 
central (r = 0) one. Indeed, a singularity cannot be naked if it occurs after the 
formation of the apparent horizon t/,(r). But using (|4.8D, we get ts{r) = th{r) 
only if \l'(r, i?/i(r)) = 0, which happens if and only if r = 0. Thus, the shell 
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labeled r > becomes trapped before becoming singular, and hence all non-central 
singularities are censored. An important consequence of this fact is that, although 
the area-radius coordinates map the singularity curve to the axis i? = 0, it is actually 
only the 'point' R = r = that has to be analyzed in order to establish the causal 
structure of the spacetime. 

Next section will be devoted to the study of the nature of the central singularity. 
A singularity is either locally naked, if it is visible to nearby observers, globally 
naked, if it is visible also to far-away observers, or censored, if it is invisible to any 
observer. We shall not be concerned here with the issue of global nakedness, so 
that our test of cosmic censorship will be performed on the strong cosmic censor- 
ship hypothesis: each singularity is invisible to any observer. We call a blackhole a 
singularity of this kind, although one could conceive situations in which the singu- 
larities are locally naked but hidden by a global event horizon, an obvious example 
of this being of course the Kerr spacetime with mass to angular momentum per 
unit mass ratio greater than one. However, in all the examples so far discovered 
of dynamical formation of naked singularities, one can easily attach smoothly to 
the region of spacetime which contains the naked singularity a regular asymptotic 
region containing no event horizon. It is, therefore, likely that the unique version 
of the cosmic censorship conjecture that can be the object of a mathematical proof 
is the strong one. 



5. The spectrum of endstates 

5.1. The radial null-geodesic equation. The equation of radial null geodesies in 
the coordinate system (r, R) is given by 

(5.1) ^ = u^/A[Y-u]. 

ar 

Indeed, all null curves in two dimensions can be reparameterized to become geo- 



desies. Therefore, ( |5.1[ ) comes from ( ^3| ) setting ds^ = 0, {d9^ + sin^ 9d(p^) = 0, 



and requiring the future-pointing character of the curve. 

The center i? = r = is (locally) naked if there exists a future pointing local 
solution Rg{r) of the geodesic equation which extends back to the singularity (i.e. 
-R(O) = 0) and "escapes from the apparent horizon", that is Rg{r) > Rh{r) for 
r > 0. We will study in full details only the existence of radial null geodesies 
emanating from the singularity. We are, however, going to prove that if a singularity 
is radially censored (that is, no radial null geodesies escape), then it is censored (see 
subsection |5]^ below). 

In what follows we shall need to consider sub and super solutions of the equation 
(|5.ip. We recall that a function yQ{r) is said to be a subsolution (respectively super- 
solution) of an ordinary differential equation of the kind y' = f{r, y) if it satisfies 
y'o < f{r,yo) (respectively >). 

5.2. The main theorem. It is easy to check that in a collapsing ARS spacetime 
(see Definition Ob the Taylor expansion of VA{r, 0) at the center is given by 



(5.2) VA(r,0) = er""' + 
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where ^ is a positive number and n is a positive integer. The following theorem 
shows that the causal nature of such spacetimes is fully governed by these two 
quantities. 

Theorem 5.1. In a collapsing ARS spacetime, the singularity forming at the center 
is locally naked ifn= l,ifn = 2,orifn = 3 and - > where 

(5.3) ^ 

Otherwise the singularity is covered. 

For sake of clarity, we divide the proof of the main theorem into two parts. 

5.3. Sufficient conditions for existence. 

Theorem 5.2. In a collapsing ARS spacetime, the singularity forming at the center 
is locally naked ifn<3orifn = 3 and ^ > ^c- 

To prove the theorem we first need the following crucial result: 

Lemma 5.3. In a collapsing ARS spacetime the apparent horizon is a supersolution 
of equation ( ]5.1D . 

Proof. This result can actually be proved for much more general spacetimes (es- 
sentially, only the weak energy condition is needed) [^. However, we give here 
a very simple proof for the model at hand. As we have seen, we always suppose 
that a is positive (see Definition ^3] ). This implies 2\['(r, 0) = ar^, and from the 
definition of the apparent horizon Rh{r), it is 

(5.4) Rhir) = 2^(r, Rhir)) = 2*(r, 0) + i?;,(r)^,^(r, Sr) = ar\ 

{Sr G (0, Rh{r))) since \I/,R(r, ^r) is infinitesimal. This implies -R^(r) = 3ar^ > 
in an open right neighborhood of 0, whereas it is easily seen that M(r, Rh{r)) = 
Y(r, Rh{r)), implying that the right hand side of ( p?T| ) gives zero when evaluated 



at (r, Rh). Hence Rh{r) is a supersolution of (^7T|). □ 



Lemma 5.4. The singularity forming at the center is naked if there exists a subso- 
lution of equation ( |5.1[ ) of the form R{r) = xr^, where x > a. 

Proof. The singularity is naked if there exists a geodesic Rg{r) such that Rg{r) > 
Rh{r) in an open right neighborhood of 0. If R{r) = xr^ is a subsolution of 
( pTTt ) and X > a, then R{r) > Rh{r). Let > and Rg{r) a geodesic through 
-Rg(ro) E]Rh{rQ),R{rQ)[. This curve cannot cross the subsolution from below, 
neither it can cross the supersolution from above. Thus Rg{r) is defined in ]0, tq] 
and \im.r^Q+ Rg{r) = 0. Therefore Rg is the sought geodesic. □ 

Theorem [5!2i We will derive sufficient conditions for R(r) = xr^ to be a subso- 



lution of (pJJ) with X > a; the result will then follow from Lemma 5.4. We have 
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(see and 
(5.6) 

since Y{0, 0) = 1. 
Now, from (^), 

(5.7) 
But 



y(r, R{r))-u{r, R{r)) = 1 



A(r,i?(r)) = VA(r, 0) - 



1 dn 



Q F(r, cr) dr 



(r, cr) dcr. 



1 dH 



r, a) da 



g F(r, cr) (9r 



° 2F(r,a)[E.+,=3V^^^ + --- 



3/2 



da 



3/2 



r dr 



r 



3/2|T=x 
lr=0 



— xr^ 



where the variable change a = r^r has been performed, and we recall that a = h^Q. 
Then 



(5.8) 



A(r, R{r)) ^ ^ r""^ + ^/ -xr" 



The curve R is certainly a subsolution of ( |5.ip if the following inequality holds: 

/a \ / fa 



(5.9) 



3xr^ < 1 



— ^ r 

a; 



+ xr^ 



This inequality holds always if n = 1 or if n = 2; namely, the term on the right 
hand side is a positive function that behaves like r™ whenever x > a. If n = 3 
condition ( Pl^j ) is equivalent to 



(5.10) 
where 
(5.11) 



S{x, |) < 0, 



S{x, p) = 2x'^ + - 3p Vtt a;^/^ + 3p a, 



and using standard techniques it can be seen that ( |5.10D holds for some x > a if 
and only if 



(5.12) 



^ 26 + 15^3 _ 
a 2 



□ 
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5.4. Necessary conditions for existence. 

Theorem 5.5. In an collapsing ARS spacetime, if the singularity forming at the 
center is locally naked then n < 3 or, n = 3 and ^ > ^c- 

To show the theorem, essentially adapting an argument exploited in for dust 
solutions, we need the following Lemma. 

Lemma 5.6. In a collapsing ARS spacetime, if a curve R{r) is a geodesic emanat- 
ing from the central singularity such that 



^_r(^) '■— / e''^H{r, a) da 
verifies \imr^Q+ t^(r) = to = limr^o+ ts(r), then 
(5.13) lim ^ = 0. 



r^0+ r 

Proof. We have \imr^Q+{ts{r) — t^{ry) = 0. Recalling that 



ts{r) = / e-'"Hir,a)da, 



we also have 



ts{r) - t^(r) 



rRir) I- 

/ e-''H(r,a)da ^ / 
Jo Jo 



B.{r)/r 



1/2 



dr, 



and since y is bounded this quantity must be infinitesimal, so that ( |5.13D must 
hold. □ 

Theorem [53[ Let R{r) = x{r)r^ be a geodesic such that x(r) > a, and -R(O) = 0. 
Using Lemma x{r)r'^ must be infinitesimal, so it is a straightforward calcula- 
tion to verify that 



(5.14) 
(5.15) 



u{r, -R(r)) 



a 



x{r) ' 

Y{r, R{r)) - u{r, R{r)) = 1 



a 
x{r) 



whereas, using the same arguments as in (|5.8[), it is 



(5.16) y/A{r,R{r)) = y/A{r,0) - [ 

Jo 



Since R{r) is a geodesic, (|5]T[) yields 
(5.17) 



r, a) da 



Y{r, a) dr 



x{r) 



-x[r)r 



a 



x'(r)r = 



a 



x(r) 



a 



x{r) 



^r"-3 + x{r) 



3 x(r). 
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By contradiction, let us first assume n > 3. Therefore, if x(r) went to +00 for 
r ^ 0+, using ( |5.17p it would be 

2 

x'(r) = — x(r)ilj(r) + (f(r), 
r 

where ip{r) is a bounded function and V^(0) = 1. Using comparison theorems for 
ODE's there would exist A > such that x(r) > ^ and so R{r) > \r, which is 



in contradiction with Lemma If limr^o+ ^i^) = h with / G (0, +00), then a 
straightforward calculation of the limit of both sides in ( |5.17D yields a contradiction 
as well, because x'{r) = —tl^ for some positive constant k. If x{r) does not have 
limit, there exists a sequence of local minima for x(r) such that ^ for 
m — > 00, x'{rm) = and x{rm) is bounded. Evaluating ( |5.17| ) for r = r^, and 
taking the limit of both sides for m — > 00 we get a contradiction once again. Then 
n < 3. In the case n = 3, equation ( |5.17| ) can be written as 

, Sixir), I) 

(5.18) x\r)r = ^ /.'^^ +..., 

x[r) 

where S{x,p) was defined in ( |5.1 1| ). Arguing as before, x{r) cannot go to +00 as 
r goes to since R{r) would be bounded from below by Ar for some A > 0. If 
limr_o+ x(r) = / < +00 then S{x{r), |) S{1, |) and, from (|TTE[), S{1, |) = 0. 



Since / > a, equation ( |5.12[ ) implies ^ > ^c- If x{r) does not have limit, as before 
there exists a sequence ^ of local minima for x{r) such that x'{r„i) = and 
x{rm) is bounded, and this yields, using ( |5.18| ), that S{x{rm), |) — > 0, that is x{r. 



converges to a root of S{x, |) = greater that a, which is possible again only if 

a — 

5.5. Non-radial geodesies. We have limited our analysis to radial null geodesies. 
However, we will now show that, if no radial null geodesic escapes from the singu- 
larity, then no null geodesic escapes at all. In other words, we show that a radially 
censored singularity is censored (in the case of dust spacetimes, this result was first 
given by Nolan and Mena [0). 

Let by contradiction R(r) be a non radial null geodesies escaping from the center, 
that we will suppose radially censored. Arguing in a similar way as for recovering 
(|5.ip, we have that R(r) solves equation 

(5.19) ^ = -u^ 

ar 



B + ^/A + CLyR^ , 

where is the conserved angular momentum. Since C is positive we have 

d^ 



dr ^ 



that is R{r) is a subsolution of the null radial geodesic equation (|5.1D. By hypoth- 
esis R{r) > Rh{r) for r > 0, and -R(O) = Rh{^)- Then a comparison argument 
in ODE similar to the one exploited in Lemma ensures the existence of a radial 
null geodesic, which is a contradiction. Thus, the following holds true: 

Proposition 5.7. Any singularity radially censored is censored. 
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We emphasize that this is a general result, not depending on the class of solutions 
we are dealing with, but only on the spherical symmetry of the model. 

5.6. Physical interpretation of the results. As we have seen, although the math- 
ematical structure of the solutions of the Einstein field equations and the way in 
which this structure governs the properties of the differential equation of radial 
geodesies is extremely intricate, our final results are nevertheless extremely simple: 
what governs the whole machinery is just the first term of the Taylor expansion of 
\/A(r, 0) near the center (Theorem |5TT| ). To understand what this results says in 
physical term, it is convenient to write VA{r, 0) as 

(5.20) VA(r,0) = /i(r)+/2(r) 

where (see also (p.9[)) 
(5.21) 



/i(r) : = 



Y{r, r) 
and we have defined 



d 

— I H{r,a)da 



7(r,a) : = 



1 



h(r) :-- 



1 



lir,cr) — {r,a) da, 



Y[r,a) Y{r,r) 
The function Ji (r) has the following behavior: 

d 



(5.22) /i(r) 



Y(r, r) dr 



i 'H{r,rT)r dr 
Jo J 



2 [P(r) - rPg(r) + . . .f^ 
where the following Taylor expansion through (0, 0) has been introduced: 

(5.23) 2^(r, R) + R{Y'^{r, R) - 1) = ^ hijr^R^ + ^ hijr^R^ + . 

i+j=3 i+j=3+p 

and 

Jo 2P(r)3/2 

The remaining summand I2 (r) is zero if Y depends only on r , i.e. if the acceleration 
A(r, R) = R) (formula ( p.22[ )) vanishes. Its behavior is as follows: 



(5.24) P(r) = J2 K-r^ Q{r 

i+j=3 



dr. 



hir) = - 



7(r,(T)- 



E 



3/2 



da 



;1 E^+j=3^^^J^' + ■ 







2r [P(r) + . . . 



dr. 



Now, using Taylor expansion of Y through the center 



F(r, R) = ifir) + ^ h^r^W + ..., 



i+j=q+l 
j>0 
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(f{r) contains the terms where R does not appear, and (^(0) = 1), it is 

'y{r,rT) ^Y{r,r)-Y{r,rT) =r^+'^ ^ kij{l - r^) + . . . = S{t) r^+'^ + . . . , 

i+j=q+l 
j>0 

where q is easily seen to be the order of the first non vanishing term of the expansion 
of the acceleration near the center. Thus, we obtain: 



• 1 

/2(r) = -br'^ + . . . , where b := I S{t) 



2P(r)3/2 



^dr. 



It follows that 

(5.25) VA(r, 0)=pa r^'^ -br'^ + ... 

The index n defined in formula is thus the smaller]] between p and g + 1. The 
value of p is clearly related to the degree of inhomogeneity of the system, since one 
can generate a low value of it taking a low order of non- vanishing derivatives of the 
initial density profile at the center (formula ( |3.17D ). This effect can be related to the 



shear as well (see e.g. [[T], [T^). The value of q, and thus the second term, is related 
strictly to the acceleration: it vanishes if the acceleration vanishes, and in any case it 
does not contributes to the nature of the final state if the system "does not accelerate 
enough" near the center. The effect of this term on naked singularities formation 
can be considered as a new feature of our models. It can, in fact, be shown that 
(virtually) all the particular cases already known in literature of formation of naked 
singularities in the gravitational collapse of continua (for instance the dust and the 
tangential stress model) can be retrieved as particular cases of our main theorem, 
and that in all such cases the acceleration term is negligible and the effect does not 



occur (of course, these cases to not exhaust the content of Theorem p7\\ ). 

6. Discussion and conclusions 

From the very beginning, the problem of Cosmic Censorship was shown to be 
linked to specific mathematical (as opposed to physically transparent) properties of 
the available arbitrary functions, like e.g. 'first derivative of initial density equal to 
zero, second non-zero at the center' . In all the cases which have been discovered 
so far the situation became more and more intricate. In the present paper, we have 
constructed a new class of solutions which contains as subcases all the solutions for 
which censorship has already been investigated in full details. The new solutions 
add to the (rather scarce) set of spherically symmetric spacetimes whose kinemati- 
cal properties are generic. The analysis of the structure of their singularities allowed 
us to show that a general and simple pattern actually exists. This pattern follows 
from Theorem ^7T]: given any set of regular data, and any equation of state within 
the considered class, the formation of naked singularities or blackholes depends on 
a sort of selection mechanism. The mechanism works as follows: it extrapolates 
the value of an integer n and selects the final state according to it. The extrapolation 
essentially depends on the kinematical invariants of the motion. If the resulting n 
equals one or two or it is greater than three, the final state is decided and has no 



'One can conceive very special cases in which the terms exactly balance each other and the index 
n has to be defined at the next order. 
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Other dependence on the data or on the matter properties. Dimensional quantities 
such as e.g. value of the derivative of the density at the center, cosmological con- 
stant, value of the derivative of the velocity profile at the center, profile of the state 
equation for tangential stresses, and so on, play a further role only at the transition 
between the two endstates, occurring at n = 3. This role is to combine themselves 
to produce a non-dimensional quantity which acts as a critical parameter. 

The (spherically symmetric) cosmic censor seems to answer to the court, that the 
formation of naked singularities or blackholes is essentially a local and kinematical 
phenomenon: it neither depends (or weakly depends) on what is collapsing, nor it 
depends on the details of the data characterizing how it starts collapsing. The 
formation or whatsoever of blackholes or visible singularities depends only on the 
kinematical properties of the motion near the center of symmetry of the system. 

Acknowledgement. The authors wish to thank the referees for their very useful 
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